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Abstract 

Collisions of heavy nuclei at very high energies offer the exciting possibility of experimentally 
exploring the phase transformation from hadronic to partonic degrees of freedom which is pre- 
dicted to occur at several times normal nuclear density and/or for temperatures in excess of 
~ 170 MeV. Such a state, often referred to as a quark-gluon plasma, is thought to have been the 
dominant form of matter in the universe in the first few microseconds after the Big Bang. Data 
from the first five years of heavy ion collisions of Brookhaven National Laboratory's Relativistic 
Heavy Ion Collider (RHIC) clearly demonstrate that these very high temperatures and densities 
have been achieved. While there are strong suggestions of the role of quark degrees of freedom 
in determining the final-state distributions of the produced matter, there is also compelling evi- 
dence that the matter does not behave as a quasi-ideal state of free quarks and gluons. Rather, 
its behavior is that of a dense fluid with very low kinematic viscosity exhibiting strong hydro- 
dynamic flow and nearly complete absorption of high momentum probes. The current status of 
the RHIC experimental studies is presented, with a special emphasis on the fluid properties of 
the created matter, which may in fact be the most perfect fluid ever studied in the laboratory. 

Key words: perfect fluid, perfect liquid, quark-gluon plasma, QGP, RHIC, ultra-relativistic heavy ion 
collisions 

PACS: 01.30.Cc 04.70.Dy 05.20.Jj 11.15.-q 11.25.Tq 12.38.Mh 25.75.-q 25.75.Bh 25.75.Ld 66.20.-d 



Preface: This manuscript is respectfully dedicated to the chair of INPC07, Prof. Shoji 
Nagamiya, for his extraordinary contributions to all of nuclear physics, but especially for 
his essential roles at RHIC, on PHENIX, at Columbia and in influencing this author's 

professional career. 

Email address: zajcfinevis. columbia.edu (W.A. Zajc). 



Preprint submitted to Elsevier 



27 September 2007 



1. Introduction 



Experiments at Brookhaven National Laboratory's Relativistic Heavy Ion Collider 
have achieved their goals of creating and characterizing a new state of matter, which 
has come to be known as the strongly-coupled Quark-Gluon Plasma (sQGP). The strik- 
ing discoveries and their imphcations from the initial three years of RHIC operations 
were extensively detailed in "white papers" from the four experiments- BRAHMS [T], 
PHENIX[2J, PHOBOSS] and STAR|5 which, together with understanding developed in 
previous theoretical work[5l, led to the announcement [B] of the "perfect liquid" behavior 
of the matter produced at RHIC. 

Since that time, analysis of substantially larger data sets has provided additional ex- 
perimental evidence in support of those statements. Perhaps more importantly, further 
discoveries have produced both new insights and new puzzles which demand more de- 
tailed experimental and theoretical investigation. Following a brief review, a sampling of 
these new experimental results and their implications will be presented. Due to length 
restrictions, the topics covered, while adhering closely to those presented at the confer- 
ence, will be even more selective than the author's talk[7j, to which the reader is referred 
for supporting material. 



2. The Initial RHIC Discoveries 



During the initial phase of RHIC operations it was quickly established that the rel- 
ative abundances and spectra of the particles produced in Au-|-Au collisions at RHIC 
energies (initially ^/sNN = 130 GeV in RHIC Run-1, followed by 200 GeV collisions 
in RHIC Run-2 in 2001-2) were consistent with emission from a thermally equilibrated 
source |:8_9_ 10, llj with a chemical freeze-out temperature Tq ~ 170 MeV and a low baryon 
chemical potential in a manner consistent with trends seen in lower energy collisions [T^]. 
However, analysis of the yields as a function of the angle cj) with respect to the reac- 
tion plane of the collision revealed the presence of strong "elliptic" flow [13^ (Figure [ij, 
parameterized in terms of the Fourier coefficient V2 in the expansion 

— - 1-f 2v2(pt)cos20+--- (1) 

(here f 2 is also regarded as a function of the transverse momentum pT = \p\ sin of 
the emitted particle). In contrast to behavior at lower energies, the strength of this 
angular modulation, and its systematic variation with the mass of the produced particles, 
was found for the first time to be consistent with the solutions of ideal (non-viscous) 
hydrodynamics. We will return to the importance of this observation in Section [3.4[ 

The discovery of "jet quenching" at RHIC|T3]- a strong suppression in the production 
of high pt particles in nuclear collisions relative to the expected yield based on p+p 
collisions, also stands in stark contrast to results at lower energy, where particle produc- 
tion at high transverse momentum is enhanced rather than suppressed in nucleus-nucleus 
collisions (Figure [iJ. Expressed in terms of the ratio Raa{pt), defined as 

^ , Yield in Au + Au events 
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Fig. 1. The discovery in y/sjvjv = 130 GeV j c Au+Au collisions at RHIC of strong elliptic flow (left, 1131 ) 
and of jet quenching (right, |14p. Left: The flow strength parameter V2 versus transverse momentum pT 
for charged particles produced at mid-rapidity in minimum bias collisions. Right: The suppression factor 
Kaa versus for """'s (circles) and charged particles (squares) in central collisions, compared to lower 
energy results. 

where the denominator consists of the p+p yield scaled, as per perturbative QCD (pQCD) 
by the equivalent parton+parton flux from a Au+Au collision, the suppression was found 
to be as large as a factor of 5 in the most central events at y^slvw — 200 GeV |15irB] . In a 
curious inversion, the realization 17J that detailed information on the opacity and other 
properties of a dense thermal QCD system could be obtained using the very deviations 
from pQCD expectations absent interactions in a produced medium spurred development 
and application of a sophisticated technologv |18ll9l20l21l22l23j making possible "tomo- 
graphic" studies of the produced matter. The observed quenching was consistent with 
parton energy loss rates ~ 15 times higher than in cold nuclear matter [27]. and demanded 
an initial matter density of order 100 times that of normal nuclear matter P4l25l26j . A 
striking observation in support of these estimates was the disappearance of the "away- 
side" jet partner in Au+Au collisions 29J (Figure [2]), indicating that the matter density 
was essentially opaque to high-pj^ partons and that the observed high transverse momen- 
tum "trigger" particles were dominated by surface emission. 

Three other early key developments can only be briefly mentioned here: 
- The interpretation of the jet-quenching results was bolstered by reliance on in situ 
measurement of baseline (p+p) and control (d+ Au) data. Comparison of the p+p data 
to theoretical calculations established the quantitative reliability of pQCD calculations 
at RHIC energies [30 . The demonstration that suppression effects were absent in d+Au 
collisions pi .32 3 3 34 provided crucial evidence that the quenching observed in Au+Au 
collisions was due to parton propagation in a dense thermal environment, rather than 
to modifications of the nuclear wave function. 
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Fig. 2. Away-side jet disappearance: the angular distribution of high pj^ (> 2 GeV/c) particles associated 
with a trigger particle p^'^ S [4,6] GeV/c measured in p+p (open circles) and Au+Au (filled circles) 
collisions as a function of centrality (expressed as the percentage of the total cross section). Taken from 
Ref. [29]. 

- While both thermal and hydrodynamic models worked with unprecedented accuracy 
to describe the final state distribution of particles at RHIC, the actual abundances 
were low compared to pre- RHIC expectations [35 . Both the low value of the charged- 
particle density [36] and its systematic variation with centrahty and energy [37)40|38j 
were described quantitatively in a model incorporating gluon saturation in the initial 
nuclear wave-function j39j. Subsequent developments have shown the importance of this 
observation for establishing reliable initial conditions for hydrodynamic calculations at 
RHIC energies m]. 

- The surprising observation that baryon and anti-baryon yields were comparable to 
those of mesons at intermediate (4-6 GeV/c) transverse momenta|12] was not under- 
standable in terms of standard thermal production models with a chemical freeze-out 
condition. However, this so-called "baryon anomaly" found a natural explanation in 
recombination models [I3l44l45l46l47j . which in a sufficiently dense partonic medium 
will favor the creation of final-state hadrons via coalescence of lower-momentum quarks 
drawn from a thermal spectrum over the fragmentation of a higher-momentum parton 
to a hadron with lower momentum. The successes of this mechanism, which clearly 
requires partonic degrees of freedom in a dense medium, led the authors of Ref. [44^ to 
note "... our scenario requires the assumption of a thermalized partonic phase char- 
acterized by an exponential momentum spectrum. Such a phase may be appropriately 
called a quark-gluon plasma." 
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3. More Recent Developments 

The steady improvement in RHIC luminosity^ and in experimental data-taking ca- 
pabilities has led to an increase approaching three orders of magnitude in integrated 
luminosity over that of the initial (Run-1) discovery period. It is not possible to enumer- 
ate, much less fully explore, all of the results obtained from this cornucopia of data. (See 
Ref. [IS] for a recent review.) Instead, focus will be applied to a very hmited subset of 
topics most relevant to future investigations that seek to characterize the properties of 
the medium. 

3.1. Direct Photon Measurements 

While the absence of suppression effects observed in d-|-Au collisions also suggested that 
the pQCD scahng methodology was well-controlled, it obviously did not directly demon- 
strate that this was also the case for Au-|-Au collisions. High energy {px >~ 3 GeV/c) 
direct photons, while experimentally very difficult to separate from the copious back- 
ground from tt" and rj decays, do provide the desired calibration, since their production 
rate is proportional to the initial parton-|-parton flux and they interact only very weakly 
with medium. The measurement of Raa for direct photons in Au-|-Au collisions [49], 
which requires tight control of experimental systematics over several orders of magni- 
tude, clearly establishes the validity of the assumed scaling techniques. 

The effort now turns to ever-increasing precision in measuring and normalizing the 
observed photon yield, in search of predicted ~ 20% effects due to the interplay of 
isospin, fragmentation, shadowing and energy loss |50|5"T] . An important experimental 
development is the first proof-of-principle results on gamma-hadron correlations |52|5"3] . 
which are a necessary precursor to the long-desired goal of using direct photons as a 
calibrated tag to measure precisely jet energy loss in nuclear collisions |54j. 

3.2. Detailed Investigations of Hydrodynamic Behavior 

The consistency of the RHIC experimental data on elliptic flow with calculations based 
on ideal hydrodynamics, together with the jet quenching results demonstrating the ex- 
traordinary density of the matter, resulted in the descriptor "perfect liquid", in analogy 
with usage of "perfect fluid" in general relativity to denote a fluid that is completely 
isotropic to co-moving observers [55 j (thereby implying zero viscosity and perfect heat 
conductivity) . A great deal of experimental and theoretical work is underway to deter- 
mine the kinematic regime in which this description (approximately) applies, and to 
quantify the transport properties of the near-perfect medium. A highly restricted sample 
of these efforts is presented in this section. 

3.2.1. Scaling Behavior of Elliptic Flow 

As noted in Section [2] data from RHIC show that the detailed dependence of the ellip- 
tic flow parameter V2{pt) on particle mass is consistent with calculations based on ideal 

^ RHIC now routinely operates at more than 4 times its design luminosity of 2 X lO'^'^cm^^s^^ for 
Au+Au collisions. 
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Fig. 3. The scaling|57] with transverse kinetic energy KEt of PHENIX and STAR l58l59l data for the 
elliptic flow parameter V2{pt) for various particles. 

hydrodynamics. This observation has been considerably sharpened by the discovery of 
a family of exact scaling solutions to ideal hydrodynamics [S5] which show a universal 
scaling behavior of the elliptic flow parameter V2 upon a reduced kinematic variable. 
Near mid-rapidity, this variable reduces (approximately) to "transverse kinematic en- 
ergy" KEx = ^y m? + — m. Indeed, the mass-dependent structure in V2{pt) does fall 
on such a scaling curve [FT, as shown in Figure [s] While the scaling behavior appears 
limited to KEt < 1-0 GeV, note that a) more than 98% of the produced particles 
are in this domain and b) the larger values of KEt correspond to a regime where hard 
scattering may already play a role, potentially explaining the observed deviations from 



hydrodynamic behavior (however, on this point see further discussion in Sections 3.4 and 
3.5). With these mild caveats, the scaling behavior shown in Figure [3] is consistent with 



a key prediction of ideal (inviscid) hydrodynamics for bulk particle production at RHIC. 

3.3. Bounding Perfection 

Taken at face value, the good agreement between the experimental data and inviscid 
hydrodynamics supports Landau's observation [BD] that the very conditions for the apph- 
cability of hydrodynamics to nuclear systems (a mean free path much smaller than the 
system size) necessarily lead to a negligible viscosity. One of the most fascinating devel- 
opments in recent years has been the conjecture [61] that there may be a fundamental 
bound from below on the value of viscosity, that is, for any thermal fluid, the uncertainty 
principle requires a non-zero viscosity. More precisely, the conjecture states that the ratio 
of viscosity ij to entropy density s must satisfy 

■ (3) 

s 47r ^ ' 
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Particularly intriguing is the origin of the bound. While the simple dependence on 
h implies it is a strictly quantum mechanical result, and while estimates based on un- 
certainty relations have been made [62i,63j . the first explicit derivation of a numerical 
value for the bound was obtained for a maximally supersymmetric Yang-Mills theory 
via the AdS/CFT correspondence [64J . By exploiting a duality between (string) quantum 
gravity in a higher-dimensional Anti de Sitter(AdS) spacetime and conformal field the- 
ory (CFT) on the boundary of AdS, strongly coupled problems in the field theory are 
mapped onto weakly coupled, and thus semi-classical, gravity calculations in the bulk. As 
a result, the entropy density in the field theory is dual to the entropy of extended black 
branes, and the viscosity is dual to graviton absorption by the brane. While the super- 
symmetric conformal field theory would appear to be far removed from non-conformal 
non-supersymmetric ordinary QCD, arguments have been made that for thermal QCD 
systems somewhat above the critical temperature may at least qualitatively be regarded 
as a conformal theory in which the only dimensionful quantity is the temperature |66|67] . 

3.4. Approach to Perfection 

Given the (conjectured) existence of a bound on ry/s, it is only natural to ask how 
closely does the QGP fluid produced at RHIC approach the bound. The most natu- 
ral approach is not to separately extract 77 and s from the data, but instead to note 
that damping (of sound, flow, etc.), at temperature T and zero chemical potential is 
proportional to This was the basis for a first schematic calculation [68, that was 
instrumental in developing the case for ideal fluid behavior at RHIC [55] . Since then, 
somewhat more sophisticated analyses coupled with greatly improved data sets have led 
to first attempts to determine the range of allowed values for 77/s in ^snn — 200 GeV 
Au-|-Au collisions at RHIC. Estimates based on the magnitude of elliptic flow indicate 
7?/s = (1.1 ± 0.2 ± 1.2)j^[7D; and r]/s = (1.9 - 2.5)^ 71J , an analysis of the damping 
of pt fluctuations gives 77/s = (1.0 — 3.8)^ [72; , while a simultaneous analysis of the 
observed energy loss and flow of heavy quarks produces rj/s = (1.3 — 2.0)^[73 . These 
different methods provide consistent support for the conclusion that the produced mat- 
ter is within a factor of 2-3 of the conjectured bound, that is, has a viscosity to entropy 
density ratio lower than any other known fluid. The above values are comparable to, 
although generally smaller than, the value 77/s ~ 0.5 ~ 6 (^) obtained from a study of 
the breathing modes of a gas of cold trapped ^Li atoms at the unitary limit [74 , which ex- 
hibits a viscosity to entropy density ratio even lower than liquid helium near the lambda 
point [61] ■ 

3.5. Reaction of the Medium 

While the jet energy loss phenomena would appear to be strictly in the domain of 
perturbative QCD, and therefore seemingly divorced from the bulk flow behavior at 
low transverse momenta, a fascinating connection between these effects has emerged 
in recent studies [75] . The disappearance of away-side jets shown in Figure [2] strongly 
suggests that the directed energy in a high momentum transfer parton scatter is absorbed 
by the medium. A sufficiently dense and strongly coupled medium, while thermalizing 
that energy, must nonetheless conserve momentum. One method for doing so is via the 
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Fig. 4. A polar representation of the data from Ref. |78l . showing the development of strongly modified 
away-side jet response as a function of centrality. 

development of a shock or "Mach" cone |76l77] . i.e., the transfer of the directed energy of 
an away-side jet into a collective excitation of lower momentum particles in the medium. 
Indeed, after carefully taking into account the modulation of the background due to 
elliptic flow, both published data[78 79J and more recent preliminary results [80 81'82j 
provide support for this hypothesis from an examination of the angular correlations of 
away-side particles in the range ^ 1 GeV/c < px <^ 2 GeV/c. 

Further work, both experimental and theoretical, is required to determine if the clearly 
observed away-side jet distortions at low px are in fact Mach cones. Experimentally, 
more detailed analysis of both 2 and 3-particle correlations in more extensive data sets is 
underway. Theoretically, moving from a bulk description of hydrodynamic shock fronts 
to a more microscopic view of the energy-momentum transport is very challenging. Here 
too AdS/CFT methods have played a role, as the duality permits calculations at all 
length scales. While this approach suffers not only from the standard concern that the 
gauge theory studied in not QCD but also from the restriction to calculating wakes from 
infinitely massive quarks, it has nonetheless provided substantial insight into the energy 
flow and medium response [83] . 

3.6. Puzzles from Quark Number Scaling 

The recent spectacular advances in our understanding of hot QCD have been driven 
by experimental discoveries, together with the subsequent development of detailed the- 
oretical understanding. It is therefore appropriate to discuss briefly at least one area 
where such understanding remains elusive. A case in point is provided by the observed 
quark number scaling of elliptic flow |57l8"l] . Examination of the two branches of ■i;2(KET) 
displayed in Figure [3] shows that the upper branch consists of baryons; the lower branch 
mesons. Upon further scaling of both KEt and t;2(KET) by Uq, the number of constituent 
quarks, the two branches merge into a universal scaling curve for elliptic flow (Figure |5|. 
A critical test of this observation is the (f> meson- although more massive than a nu- 
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Fig. 5. The additional scaling of V2(KEj') with the number of constituent quarks n„|57|. 

cleon, it follows the scaling curve for mesons, not baryons[85 86 . While it is tempting to 
conclude that the strict scaling according to constituent quark content provides incon- 
trovertible evidence for the underlying role of quark degrees of freedom in establishing 
the elliptic flow, such a conclusion appears to be at odds with the observation of perfect 
fluidity. Quasi-particles with lifetimes comparable to the relatively long time over which 
hydrodynamic flow persists are incompatible with the very short mean free paths (and 
hence large widths and short lifetimes) implied by low viscosity f 87'. This fact, already 
noted in the earliest calculation of viscosity via the AdS/CFT correspondence [55] . has 
yet to be reconciled with the strong empirical evidence for quark number scaling of flow 
phenomena. Recent calculations [88^ which indicate that the viscosity at RHIC may be 
even lower than the conjectured bound only sharpen the dichotomy. 

4. Outlook 

The RHIC experiments have conclusively demonstrated the feasibility, and indeed the 
desirability, of performing sensitive measurements of heavy ion phenomena in a collider 
environment. The discoveries made in the initial years of RHIC operation have initiated 
a paradigm shift in our understanding of deconflned QCD matter in the region of the 
transition temperature. In doing so they have also set the stage for the heavy ion program 
at the LHC, where Pb+Pb collisions at y^'sNN = 5.5 TeV are planned following heavy 
ion commissioning in 2009. This unprecedented increase in center-of-mass energy will un- 
doubtedly result in a new round of discoveries, to be made by the ALICE experiment |89]. 
which is dedicated to the heavy ion program, together with the two large p-l-p collider 
detectors ATLAS[91 and CMS^90j. Questions of immediate interest are the fundamental 
properties of the QGP at the expected higher initial temperatures at the LHC- Will it 
remain strongly coupled? Will it exhibit ~perfect flow? However, the RHIC experience 
has taught us that the very nature of the important questions changes in response to 
experimental data, and it would be wise to anticipate that the same will be true for 
heavy ion physics at LHC energies. 
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A compelling program of more detailed investigations at RHIC is already underway, 
based on a series of detector upgrades to PHENIX and STAR, together with ongoing 
and future substantial increases in RHIC luminosity. The unique properties of bulk, 
thermalized QCD matter at RHIC energies will be better quantified via increasingly 
differential measurements, such as associated particle production in direct photon events, 
jet tomography with respect to the reaction plane, and the suppression and flow patterns 
of both charmonium and open heavy flavor. 

A particular strength of the RHIC facility is the wide range of collision energies, from 
the top energy of ^/snn ~ 200 GeV to values as low as a few GeV. This, together 
with the collider geometry which provides excellent control of acceptance systematics 
as the center-of-mass energy is varied, have stimulated interest in the search for the 
QCD critical point j92j at RHIC jS^ . Especially intriguing in this regard is the possible 
connection |94j between physics in the region of the critical point and the minimal value 
of the viscosity to entropy density ratio. The simultaneous progress on multiple energy 
frontiers at RHIC, at the LHC and at GSI/FAIR 95J are certain to lead to new insights 
and a deeper understanding of the fluid nature of quark-gluon plasma. 
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